Matrices coupled in a chain. I. Eigenvalue correlations 
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Abstract. The general correlation function for the eigenvalues of p complex hermitian 
n x n matrices coupled in a chain is given as a single determinant. For this we use a slight 
generalization of a theorem of Dyson. 
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1. Introduction 

The probability density exp[— tr V(A)] for the elements of a n x n matrix A is known to 
give rise to the probability density [1] 



F(xi, . . . , x n ) oc exp 



n 



l<i<j<n 



for its eigenvalues xi, a; n . Here V(x) is a real polynomial of even order, the coefficient 
of the highest power being positive and (3 is the number of real components of a general 
element of A, i.e. /3 = 1 if A is real symmetric, (3 = 2 if A is complex hermitian and /3 = 4 
if A is quaternion self-dual. 

The case of coupled matrices may be represented by a graph where each matrix is 
represented by a point, and two points representing matrices A and B are joined by a line 
if the coupling factor exp[c tr(AB)] is present in the probability density. When several 
matrices are coupled the probability density for the eigenvalues is known only in the case 
where these matrices are complex hermitian and the graph has a tree structure, i.e. does 
not have a closed path. 

In what follows we will consider the simplest case of a tree, i.e. that of a chain of p 
complex hermitian n x n matrices with the probability density for their elements 



F(A, 



A p ) oc exp 



-tr 



V 1 (A l )+V 2 (A 2 ) + --- + V p - 1 (A T 



-i) + \v p (Ap) 



x exp [tr {ciAiA 2 + c 2 A 2 A 3 H h c p _iA p _iA p }] 



(1.2) 



where Vj(x) are real polynomials of even order with positive coefficients of their highest 
powers and Cj are real constants. For each j the eigenvalues of the matrix Aj are real and 
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will be denoted by xj := {xji,Xj2, ...,Xj n }. The probability density for the eigenvalues of 
all the p matrices resulting from Eq. (1.2) is [2-5] 



where 



F(xi;...;x p ) 

= C exp 



n [I 1 ] 

^2 ) 2 Vl ( Xlr ^ + V ^( X 2r) H h Vp-xiXp-Xr) + - V p (x pr ) \ 

r=l ^ - J 



l<r<s<n 

x det [ e ClXl - X2s ] det [ e c ^ x ^] ■ . -det [ e c *>-i^-i^ 

ii-" 1 



c 



l<r<s<n 



det [w fc (x fcr , a; fc+ i a )] rja=li _ 



k=i 



Wk{£,v) : = exp 



■^(0-^+1(7/) + ^ 



(1.3) 
(1.4) 

(1.5) 



and C is a normalisation constant such that the integral of F over the np variables Xi r is 
one. We will be interested in the correlation functions 



Rki,...,k p ( x lli 2-lfcii •••> ^pl) x pkp) 
P V 1 n 



= 1 v J; r,-=fc,-+l 



(1.6) 



This is the density of ordered sets of kj eigenvalues of Aj within small intervals around 
Xji, Xjkj for j = 1, 2, p. Here and in what follows, all the integrals are taken over 
—00 to 00. 

The case p = 2 was considered earlier [6,10]; the expressions given in [6], though 
correct, can be put in a much simpler form: the general answer can be written as a single 
m x m determinant with m = k\ + + • • • + k p . The result is given in section 2 and 
the proof in section 3. Our result is a generalization of Dyson's for a single hermitian 
matrix [7] , according to which the correlation function of k eigenvalues is given by a k x k 
determinant: 

i2fc(xi,...,x fc ) = det[K(xi,Xj)] iJ=li k . (1.7) 



2. General correlation function 

To express our result we need some notations. Recall that a polynomial is called monic 
when the coefficient of the highest power is one. With a monic polynomial Pj(£) of degree 
j let us write 

Pijit) ••= Pj(t), (2-1) 
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and recursively, 

Py(£):= j Pi-iMm-i(v,Z)dri, 2<i<p. (2.2) 
Similarly, with a monic polynomial Qj(£) of degree j we will write 

Q,AO: QM). (2.3) 

Qij(0-= J ' Wi(Z,ri)Q i+lj (ri)dri, 1 < i < p - 1. (2.4) 

With arbitrary monic polynomials -P,(£) and Qj(£) of degree j, j = 0, 1, 2, we can 
write the product of differences asnxn determinants 

II & - = det K _1 ] = detft-ite)] = detlQ^!^)]. (2-5) 

l<i<j'<n 

Now if the polynomials Vj and the constants Cj are such that the moment matrix [M^-], 
i,j = 0, 1, n is positive definite for every n, where 



J C (w 1 *w 2 * ...*tw p _i) (tv) V J <%dri, (2.6) 



and 



(w^ *w i2 *-- -*Wi k )(Z,Ti) := J w il (£,£ 1 )w i2 (£ 1 ,&)--w ik (£ k - 1 ,r])d£ 1 ...d£ k - 1 , (2.7) 

then it is always possible to choose the polynomials Pj(£) and Qj(£) such that 

J Pj{g){w\ *w 2 * ... * v)Qk(v)d£di] = hj8 jk , (2.8) 

i.e. they are orthogonal with a non-local weight. This means that the functions Pij (£) and 
Qij(£), which are not necessarily polynomials, are orthogonal 

J Pi j (Z)Qik(Z)dt = h j 6 jk , (2.9) 



for % = 1, 2, p and j, fc = 0, 1, 2, 
Now define 



where 



Kij&ri) -Hij&ri-Eijfori), (2.10) 



n— 1 1 

0, if % > j, (2.12a) 

E ij (^r J ):={ Wifori), Xi = j + 1, (2.126) 

(w i *w i+1 *'-'*w j - 1 )(^,r]), ifi<j + l. (2.12c) 



Theorem. The correlation function (1.6) is equal to 

Rk\ ,...,k p (^ll j •••■> %lki i •••> ^-pfcp ) 

= det [Ky (xi r , ^js)] i)J= i ) ... ) p ;r= i ) ... )fci;s= i v .. )fcj • ( 2 -13) 

This determinant has m = k\ + • • • + k p rows and m columns; the first k\ rows and k\ 
columns are labeled by the pair of indices lr, r = 1, fci; the next k<i rows and k2 
columns are labled by the pair of indices 2r, r = 1, ki and so on. Each variable Xi r 
appears in exactly one row and one column, this row and column crossing at the main 
diagonal. If all the eigenvalues of a matrix Aj are not observed (are integrated out), then 
no row or column corresponding to them appears in Eq. (2.13). If all the kj are zero, then 
the correlation function is 1. 



3. Proofs 

The theorem is a consequence of the following two lemmas. 

Lemma 1. The np x np determinant det[Kij(xi r ,Xj 8 )], i,j = l,...,p; r, s = l,...,n, is, 
apart from a constant, equal to the probability density i^Xx; ...;x p ), Eq. (1.4), 

det[K ij (x ir ,x ja )] i j = lt ^ p = ( JJ^ 1 ]C- 1 F(x 1 ;...;x p ). (3.1) 



r, s = 1, n 



,1=0 



Lemma 2. Using the obvious notation of Eq. (2.7), 

(/*</)(£, »7):= / f(Z,Og((,v)d(, (3.2) 

let the p 2 functions Kij(x, y), i,j = 1, ...,p, be such that 

{K ik , if i > j > k, 
-Kik, if i < j <k, (3.3) 
0, otherwise. 

Then the integral of the m x m determinant det [K^ (xi r , Xj S )], = 1, ...,p; r = 1, 
s = 1, kf, ki > 0, k p > 0; m = k± + k<i + ... + k p ), over X£ t is proportional to the 
(m — 1) x (m — 1) determinant obtained from it by removing the row and the column 
containing the variable X£ t - The constant of proportionality is ot£ — k^ + 1, with 

a t = j Ku{x , X )d X . (3.4) 



Let us recall here a result of Dyson [7,8]. 
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Dyson's theorem. Let the function K(x, y) be such that K * K = K, then 



/ 



det[K{x h Xj)] iJ=1 ^_^ n dx n = (a - n + 1) det[K(xi, Xj)]i,j=i,-,»-i» 



(3.5) 



with 



/ 



a = K{x, x)dx. 



(3.6) 



Our lemma 2 above is a generalization of this one when the matrix elements K{xi, Xj) are 
replaced by ki x kj matrices Kij(xi r , Xj s ). 



Proof of lemma 1. Consider the np x np matrix [H^ (x ir , Xj 8 )], Eq. (2.11), the rows of 
which are denoted by the pair of indices ir and the columns by js; i,j = 1, p; r,s = 1, 
n. This matrix can be written as the product of two rectangular matrices [Qa( x ir)] 
and [Pje(xj S )/hi] respectively of sizes np x n and n x np, with I = 0, 1, n — 1. The 
rows of the first matrix [Qu(xi r )] are numbered by the pair ir and its columns by £. For 
[Pji(xjs)/h>e] the rows are numbered by t and the columns by js. Cutting the matrix 
[Hij(xi r , Xj S )] into n x n blocks, we can write 





[QiPi 


Q1P2 ■ 


■ Q\P P ' 


H = 


Q2P1 


Q2P2 ■ 






-QpPi 


QpP2 ■ 


QpPp - 



Qi 
Q2 

Qp 



[Pi p 2 



npXnp 



P inxnp 



(3.7) 



npXn 



where [<2i] r ^ := [Qu(xir)] an d [-Pj]^ s := [Pji( x js)/h(] are n x n matrices. The rows of 
Pi and the columns of Q p contain distinct monic polynomials, their ranks are therefore n. 
Thus the rank of [Hij (xi r , Xj S )] is n. As seen from Eq. (3.7) above, the first n columns of 
this matrix are linearly independent, while the remaining n(p — 1) columns can be linearly 
expressed in terms of the first n columns. In view of Eqs. (2.10) and (2.12a) the first 
n columns of [Hij(xi r , Xj S )] are identical with the first n columns of [Kij(xi r , Xj S )]. The 
determinant of the later is therefore not changed if we subtract from its last n(p — 1) 
columns the corresponding n(p — 1) columns of the former. Thus 



dot ^-Z^T^j i^^'if 7 



JSJ 



— det [ Hn {xi r , x\ 8 ) Eij(Xi r , Xj s )\ ^ _ ^ p. j _ 2 p- 

r,s = 1,2, ...,n; 



det 



r 

H21 

Hp-i 1 
- Hpi 



-E12 —E13 
— E23 













-E 
-E 



i P 

2p 



-E 1 
j -'p—i p 





(3.8) 
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From Eq. (2.12a) the last n rows of this matrix corresponding to % = p have non-zero 
elements only in the first n columns; also the matrix [Eij(xi r , Xj s )] is block triangular, 
Eij(£,r)) being zero for i > j. Therefore, 

det [Kij(xi r , Xj s )] = det [H p i{x pr , xi s )\ det [Eij{Xi r , Xj S )\ j = \^ ...,p — 1; j = 2, ...,p; 

r,s = 1,2, ...,n; 

v 

= det [H pl (x pr ,x ls )] JJdet [E j - lj (x J -i r , x js )] 
i=2 

'n-l \ p-1 

JJ /i^ 1 J det [Qpe(x pr )] det [Pi^(xi s )] JJ det [^-(xj,., x i+ i s )] , 

(3.9) 

and from Eqs. (1.4) and (2.5) one gets Eq. (3.1). This ends the proof. 

The learned reader will have recognized that the above E^s play the same role as the 
e in Dyson's proof in the case of a single matrix of the circular orthogonal ensemble [7,9]. 

Proof of Lemma 2. We want to integrate the mxm determinant [Kij(xi r , Xj s )] over xgt- 
We can write the expansion of the determinant as a sum over ml permutations, writing 
these permutations as a product of mutually exclusive cycles. The variable xm occurs in 
the row and the column labeled by it (recall that rows and columns are labeled by a pair 
of indices). If it forms a cycle by itself, then by Eq. (3.4) integration over xg t gives a 
factor ai, an d its coefficient is just the expansion of the (m — 1) x (m — 1) determinant 
obtained by removing the row and the column containing xgt- If it occurs in a longer 
cycle, say in the permutation a = (ir, it,js, •••)('")'"> then from Eq. (3.3) integration 
over X£ t decreases the length of the cycle containing It by one, giving the permutation 
a' = (ir,js, •••)(•••)"■) an d multiplies it by a factor +1 if i > i > j, by — 1 if i < i < j 
and by otherwise. Also, the cycle lengths differing by one, the permutation a has a sign 
opposite to that of a'. So the question is, given the permutation a', in how many ways 
can one insert it in any of its cycles with the algebraic weights +1, —1 and to get the 
permutation a. 

Let us represent the cycles of permutations by a graph. Only the first indices i, 
j, ... of each pair of indices will be plotted against the place number where they occur. 
For example, the cycle (42,26,36,15,24,22) is represented on figure 1, where points at 
successive heights 4,2,3,1,2,2 are joined successively by line segments or "sides". Note 
that we identify the 7 th point and the first one. Permutation a' is thus represented by a 
certain number of closed directed polygons corresponding to its mutually exclusive cycles. 
Addition of it in one of the cycles of a' amounts to the addition of a point at a height 
i in the corresponding polygon. If this added point lies on a non-ascending side, then 
the coefficient multiplying the corresponding a is +1, if it lies on an ascending side, the 
coefficient is —1, and this coefficient is zero otherwise. In other words, each downward 
crossing of the line at height i, with or without stops, contributes a factor +1, each point 
on this height contributes a factor +1 and each upward crossing, with or without stops, 
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contributes a factor —1. The graph of a' consisting of closed loops, the number of upward 
crossings is equal to the number of downward crossings at any height. 



4 




3 



2 



k=2 



Figure 1: the permutation a' = (4, 2, 3, 1, 2, 2). 



The algebraic sum of all such coefficients is thus seen to be the number of points at height 
I in the graph of a', i.e. it is kg — 1. Also the permutations a and a' have opposite signs, 
since only one of their cycle lengths differs by unity. Thus 



det[K ij (x ir ,x js )} mxm dx et = (a e - k t + 1) det [K^ (x ir , ^)] (m _ 1)x(m _ 1) , (3.10) 



where the integrand on the left hand side is a m x m determinant, i,j = 1, p; r = 1, 
s = 1, kj] m = ki + k<i + ■ ■ ■ + k p and the result on the right hand side is the 
(m — 1) x (m — 1) determinant obtained from the integrand by removing the row and the 
column containing the variable xm- This ends the proof. 

Using Eqs. (2.9-2.12), one verifies that the Hij and Eij satisfy the following relations 



Hij * Hjk 



(3.11) 




(3.12) 



(3.13) 



(3.14) 



This implies for the Kij the relations (3.3). Also 




(3.15) 
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Using lemma 1 once and lemma 2 several times one gets the normalization constant C, 

n-l 

C={nl)- p Y[hJ 1 . (3.16) 

£=0 

Again with repeated use of lemma 2 and from the definition, Eq. (1.6), of the correlation 
function Rk 1 ,...,k p (x 11 , xi kl ; ...; x pl , x pkp ), one gets Eq. (2.13). 



4. Conclusion 

The correlation functions of eigenvalues of a chain of random hermitian matrices can thus 
be written in a very compact form as a single determinant. This result may be used to 
study the large n limit of correlations between eigenvalues [10,11]. 
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